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Abstract 

In the past, the possibility to employ (scalar) material reference systems in order to describe 
classical and quantum gravity directly in terms of gauge invariant (Dirac) observables has been 
emphasised frequently. This idea has been picked up more recently in Loop Quantum Gravity 
(LQG) with the aim to perform a reduced phase space quantisation of the theory thus possibly 
avoiding problems with the (Dirac) operator constraint quantisation method for constrained 
system. 

In this work, we review the models that have been studied on the classical and/or the quantum 
level and parametrise the space of theories so far considered. We then describe the quantum 
theory of a model that, to the best of our knowledge, so far has only been considered classically. 
This model could arguably called the optimal one in this class of models considered as it displays 
the simplest possible true Hamiltonian while at the same time reducing all constraints of General 
Relativity. 
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1 Introduction 



In series of seminal papers [U [21 El S] j Kuchaf and his collaborators have constructed a whole 
class of manifestly spacetime diffeomorphism invariant matter actions which have the remarkable 
feature to clarify the conceptual setup of General Relativity. A further model of this type, 
originally introduced by Rovelli and Smolin in [5], was analysed by Kuchaf and collaborators 
in [6]. The motivation for introducing those matter actions is to use the matter considered as 
a material reference system. This allows to isolate the pure gauge degrees of freedom (under 
spacetime diffeomorphisms) encoded in general relativity in an elegant way and yields to theories 
where gauge invariant (Dirac) observables can be constructed by means of those matter reference 
fields. The matter considered in [U El E] was coined dust matter because its energy momentum 
tensor is that of a pressure free perfect fluid. The advantage of the availability of such a 
description of the phase space of matter and geometry with regard to quantum gravity was 
already emphasised in those works although a concrete quantum framework was not available 
at that time and the quantum theories following from this models could only be discussed at a 
formal level. 

A particular feature of the matter models considered is that they lead to a deparametrised 
form of general relativity (plus, if considered, standard model matter) and that the resulting 
hypersurface deformation algebra for the Hamiltonian constraints, that is the algebra among the 
individual constraints, becomes Abelian. This can be achieved by writing those Hamiltonian 
constraints in an equivalent form in which they are linear in the reference matter momenta. 
This is an important property because in this form the constraints form a true Lie algebra in 
contrast to the Hamiltonian constraints in its standard form. As far as the quantum theory is 
considered, this is an advantage because many techniques that aim at solving the constraints in 
the quantum theory (such as group averaging [7]) cannot be applied when there occur structure 
functions instead of structure constants. 

The equivalent form of the Hamiltonian constraint that has mutually commuting Poisson 
brackets with itself is of the typical form C tot (x) = P(x) + h(x) where P is the momentum 
conjugate to one of the dust fields T (called the clock field) and h is a scalar density of weight 
one built from the non dust, spatial scalar density weight one contribution C to the Hamiltonian 
constraint, the (square root of) the determinant Q = \J det(q) of the spatial metric q a b, the 
density two spatial scalar D = q ab C a Cb where C a is the non dust contribution to the spatial 
diffeomorphism constraint, the density one respectively zero spatial scalars V = q ab T ja Cb, U = 
q ab T a Ti,. For a certain subclass of models, the dependence of h on V, U is trivial and it then 
follows from the first class property of the constraint algebra that the h(x) themselves must be 
mutually vanishing. 

In [8] the general solution of phase space functions h(x) of density weight w built from 
C(x), Q(x), D(x) (but not on V(x), U(x)) with mutually commuting Poisson brackets was found 
in terms of a first order PDE which can be solved by the method of characteristics, thus extending 
the set of solutions that were obtained from the concrete models [H El El [6] but without 
providing concrete covariant actions from which they result. In [9] it was shown that all of these 
solutions can be obtained by using covariant actions of the Lagrange multiplier type used in 

m m ei m um- 

On the other hand, one can ask whether there are covariant scalar matter models that do 
not rely on Lagrange multipliers at all. This question was studied in |10| and it was pointed out 
that provided that the scalar field has pure derivative coupling to the gravitational field then 
such models again give rise to mutually Poisson commuting Hamiltonian densities h. Also the 
inverse question can be analysed using the methods developed in fit)], namely, given a general 
solution h as described in [8] (with density weight one), what is the Lagrangian L that gives rise 

1 I.e. they define the same constraint surface and the same Dirac observables. 
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to it? The solution, when it exists, can be written implicitly as the solution of an ODE. 

The general framework for describing the reduced phase space of a constrained system with 
a generally covariant Lagrangian in manifestly gauge invariant form has been considered, to the 
best of our knowledge, for the first time in I12j . Aspects of this have been rediscovered in 
different language by several authors, e.g. |13l [T4l 115] . In these works, a complicated formula 
appears that "projects" a non gauge invariant function / on phase space to a (weakly) gauge 
invariant one. It involves an infinite series of multiple Poisson brackets between / and the 
Hamiltonian constraint which makes it apparently impractical to use. However, as pointed out 
in [16] . the Poisson algebra of these formal observables can nevertheless still be simple, thus 
enabling in principle a concrete reduced phase space quantisation approach to quantum gravity. 

In \17\ [T8] the frameworks of [TJ [PL] [T2] were combined in order to compute the classical 
reduced phase space of pQ. Namely, in [U II HI HI H] the reduced phase spaces where only 
described for some models in suitable gauges and the last missing step before arriving at the 
fully reduced phase space was to use the projector formula provided in (111 [T2] and construct 
general observables with respect to the Hamiltonian constraint. 

The reduced phase space described in [T71 [18] was then quantised using the technology of 
Loop Quantum Gravity [19\ [20] in [21] and thus providing for the first time a concrete model 
for LQG where all constraints are solved. This means that all operators in question are Dirac 
observables and time evolution is driven by a physical (i.e. non vanishing, gauge invariant) 
Hamiltonian for which a concrete quantisation was provided. In [21] two possible quantisations 
of the reduced physical phase space were discussed. In the first one the usual kinematical Hilbert 
space of LQG becomes the physical Hilbert space, no solutions of the constraints have to be com- 
puted, the usual LQG inner product is the physical inner product here. In the second possibility 
techniques from the Algebraic Quantum Gravity framework introduced in [26] are considered 
and the quantum theory is defined on an abstract algebraic graph using von Neumann's infinite 
tensor product Hilbert space. 

Later, in [22J, a scalar matter model of the type considered in |10j . namely the model of [3 [6] was 
used in order to perform the reduced phase space quantisation by methods from LQG of General 
Relativity with a matter content different from [21]. While in [21] four scalar reference fields 
are considered the model in [22] contains only one scalar field. The most significant difference 
between these works is that the spatial diffeomorphism constraint is solved already classically in 
[2T] using three out of the four reference fields while in [22] it needs to be solved at the quantum 
level. 

A reduction of the diffeomorphism constraint at the classical level has the advantage that one 
can avoid ambiguities that can occur in the quantum reduction. These exists also for the spatial 
diffeomorphism constraint and are in general not unproblematic because those could possibly 
present an anomaly of the algebra of observables [23]. It turns out that the classical reduc- 
tion with respect to the spatial diffeomorphism constraints of the dust model of [1] has the 
consequence, when using the standard LQG representation, that the physical Hamiltonian of 
the theory (a diffeomorphism invariant quantity) needs to be implemented as an operator that 
preserves the underlying graph on which the quantum states are definecll- If one wants to avoid 
a graph-preserving quantisation which was one of the motivations for the quantum model intro- 
duced in [22], then within LQG one is forced to solve the spatial diffeomorphism constraint at 
the quantum level as shown and implemented for the first time in [25] . In this case the operators 
involved are not knot class preserving. 



2 This simply follows from the result that spatially diffeomorphism invariant operators that are graph-modifying 
are not densely defined in the standard kinematical representation of LQG as shown in in [31] 
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As far as semiclassical aspects of the models are concerned, there are two issues to be keppt 
in mind. First of all, semiclassical tools have not been developed at the level of the Hilbert 
space of spatially diffeomorphism invariant states at all, neither in the knot class preserving case 
nor in the non preserving case. Secondly, at the kinematical level, only for graph-preserving 
quantisations do exist semiclassical tools. These have been used to establish the correctness of 
the semiclassical limit of the dynamics for short time scales in [21] using techniques developed 
in |26| . However, for graph-modifying dynamical operators that occur the model in [22] one 
first needs to develop semiclassical states first for diffeomorphism invariant states to begin with 
and then in addition for knot class modifying operators before the semiclassical sector of the 
dynamics can be investigated. Notice that on the other hand in the Algebraic Quantum Gravity 
(AQG) framework of [26J here is no need to preserve algebraic subgraphs on which the algebraic 
spin network functions depend. 

Another more recent model discussed in [24] that was also already described in the appendix of 
PQ as a special case of the model in [1] results by switching off three Lagrange multiplier fields 
in the general Lagrangian of [T] which leads to a model with only one Lagrange reference matter 
field and likewise to [22] deparametrises only the Hamiltonian constraint but does not rewrite 
the spatial diffeomorphism constraints in equivalent form that they are linear in the momenta 
of the reference fields as in [21] . 

The authors nevertheless claim that they can treat their model with the methods from [26] 
without having to care about the spatial diffeomorphism constraint anymore and nevertheless 
describe the physical sector of the theory. Moreover, they claim that they have found the most 
promising classical foundation for the reduced phase space quantisation of General Relativity 
using methods from LQG in the sense that the physical Hamiltonian simplifies in the maximal 
sense. As we will see, some of these claims are too strong. 

It transpires that there are several models of General Relativity, the standard model matter 
and certain additional matter in the literature which lead to different reduced phase spaces and 
different physical Hamiltonians with different properties. It would therefore be desirable to start 
from a common platform in the form an action principle, to study the different models from a 
unified viewpoint and to analyse their different physical properties and promises. This is the 
objective of the present paper. 

This paper is architectured as follows: 

In section two we review the known results from the literature and then construct a general 
Lagrangian that encompasses and extends all models considered so far in the literature with 
minimal coupling of the metric (no derivatives) to matter. It serves as the desired common 
platform or parametrisation of theory space. Different models can be reached by switching on 
and off certain parameters. Basically the models considered so far fall into two types: I. Those 
which just deparametrise time and II. those which deparametrise all of spacetime. We describe 
their Hamiltonian analysis in a unified fashion. 

In section three we reconsider the model of type II studied in [25j and iterate arguments well 
known in the literature (e.g. [27} [28] and references therein) that reveal that some of the claims 
made in |24| are not tenable. It is is not possible to simply drop the spatial diffeomorphism 
constraint and quantising using the algebraic quantum gravity framework as claimed by the 
authors. Rather one has to treat the quantisation of the dynamics in that model more carefully. 
A possible way to deal with the spatial diffeomorphism constraint in the algebraic framework 
is to quantise it by the Master constraint method, see [27] for a discussion, if one wants to pro- 
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ceed along [26]. In |22j the reduction of the spatial diffeomorphism constraints at the quantum 
level was performed via the refined algebraic quantisation technique [7J by means of a rigging 
map yielding the usual (spatially) diffeomorphism invariant Hilbert space of LQG. The physical 
Hamiltonian was then quantised on the diffeomorphism invariant Hilbert space. 

Moreover, a careful treatment shows that there is an important issue with the choice of sign of 
a square root that has not been mentioned in [23]. Namely, if one insists on the usual energy 
conditions for the energy momentum tensor, then the sign of the physical Hamiltonian must be 
constrained in that model whence the physical Hamiltonian is still not without a square root as 
claimed, since it is constrained to equal plus or minus its absolute value. We also show that the 
model in [23] corresponds to a subsector of the model |21] which also explains why there is still 
an absolute value (and thus a square root) involved in the physical Hamiltonian of [24]. The 
fact that this is a subsector of [21] corresponding to vanishing spatial diffeomorphism constraint 
makes it also transparent that the latter cannot be dropped. 

In section four we work out the details of the LQG quantisation of the model [3j following 
the procedure of [171 IS EE] • This model could arguably called the optimal one in the class of 
models considered so far in the sense that not only also the spatial diffeomorphism constraint 
is already reduced classically, thus avoiding the ambiguities and potential problems pointed out 
in [23] in connection with the scalar product on spatially diffeomorphism invariant states, but 
furthermore, this time there is really no square root involved in the physical Hamiltonian. 
In section five we summarise and conclude. 

2 Theory Space 

There are roughly two types of scalar field models that have been studied in the literature and 
that were found convenient for investigations in quantum gravity. Models of type I are based 
on Lagrangians that involve beside the fields that play the role of the reference fields additional 
fields that couple to gravity without any derivatives. We will call the latter Lagrange multiplier 
fields. Models of type II do not use Lagrange multiplier fields but consider only reference matter 
fields in their Lagrangians. Note that both types consider only scalar fields that are minimally 
coupled to gravity and hence the model of a conformally coupled scalar field in [29] is not con- 
sidered in the classification here. 

Common to both types of Lagrangians is that the non Lagrange multiplier fields enter the La- 
grangian without a potential, that is, it depends only on the first derivatives of the fields. Those 
fields are used to construct the gauge invariant expressions of the remaining fields (geometry 
and standard matter) and as far as the final observables are constructed those fields have been 
absorbed. This can be understood likewise to the case of the Higgs mechanism of the standard 
model where three out of the four Higgs fields are 'eaten' by the boson fields yielding to a gauge 
invariant description of massive vector bosons. The main difference to the standard Higgs mech- 
anism is that here all reference field considered in the models are 'eaten' by the other degrees 
of freedom, whereas in the ordinary Higgs mechanism one scalar field remains in the theory 
and is not absorbed. The reason why no potentials of the reference fields are considered is that 
one wants the physical Hamiltonian to be independent of the additional scalar fields introduced 
because those will play the role of the reference time and spatial reference points. Hence without 
a potential, on achieves that the physical Hamiltonian is then not explicitly dependent of the 
reference time and thus the reduced physical system becomes conservative which has obvious 
advantages. In principle, however, also potential terms could be treated if one gives up on the 
requirement of a conservative system. We will review both types of models I and II below and 
point out some of their most important physical differences. 
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2.1 Type I 

The first class of models are those considered in [TJ [2j [3] and we will briefly discuss their La- 
grangians here. The famous Brown - Kuchaf timelike dust Lagrangian [lj is given, in our 
notation, by 

Ctd = ~V\ det(g)\ p [g^U^U, + 1], Up = V M T + W 3 V^S 3 (2.1) 

which depends on two pairs of respectively 4 scalar fields (p > 0,Wj) and (T,S 3 ) respectively. 
Here latin indices j, k, .. run from one 1 to 3 and greek indices p, v, ... run from to 3. We choose 
our signature convention for the spacetime metric tensor g^ u to be (—,+,+,+). The Bicak - 
Kuchaf null dust Lagrangian [2] is given by 

£ND = -^V|det( 5 )| (7^1^, U^ = W 3 V^S 3 (2.2) 

and results from (|2.1f) by setting p = 1, T = and dropping the cosmological constant term 
—p\/\ det(g)|/2. The opposite limit was taken in a discussion in the appendix of [1] where one 
sets Wj = and thus drops the dependence of (|2.ip on S 3 which results in 

£nrd = -\p\/\ det( 5 )| [g^U^U,, + 1], C/ M = V M T (2.3) 

This possibility of non rotational dust can be understood as a special case of the timelike dust 
model. 

The variables (p,T) and (Wj,S 3 ) appear symmetrically in (12. ip in the sense that if we write 
Ao := yfp, Xj := yfpWj and S° := T then (|2.ip becomes 

£td = -~V\ det(g)\ [g^U^U u + 1], U„ = A Q V M 5 a (2.4) 

The Gaussian dust model suggested in [3] treats them unsymmetrically and can be written in 
the form 

£ GD = -Ipy^detM bT (V M T) (V„T) + 1] - ^dit^ 5^(V M T)K, ^ = W^V M ^ (2.5) 

2.2 Type II 

The first type of model of this class is simply the massless Klein-Gordon action considered in [5] 
and whose canonical structure was discussed in detail in [6]. It leads, for general reasons that 
we review below, to an Abelian algebra of Hamiltonian constraints. As pointed out in [lOj . the 
mechanism responsible for this is by far not restricted to the Klein Gordon action but works for 
a general Lagrangian of the form 

c s = V\^)\L(i), j = -I^[v M r] [v,r] (2.6) 

where L is any function of the argument indicated. The first observation is that ()2.6|) in contrast 
to the type I Lagrangians depends on only a single scalar field T rather than several ones T, S J . 
The reason for this is that when one simply adds Lagrangians of the type (|2.6p or considers 
more generally Lagrangians depending on L{I\,.., Jjv) where Ik is as in (|2.6p but T replaced by 
a field Zk, then the physical Hamiltonian no longer deparametrises. It is for this reason that 
the models of type I either do not depend on the S 3 or if they do then only in the combination 
V/j, = WjV^S 3 involving the Lagrange multipliers. We will see below why this is the case. 

In any case, when we assume that observables are constructed using the reference matter 
fields and no geometric degrees of freedom we see that models of type II in contrast to some 
models of type I cannot lead to a fully reduced phase space at the classical level but only a 
reduction with respect to the Hamiltonian constraint can be performed. 
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2.3 Poisson commuting Hamiltonian constraints 

One particular feature of the models in [U [H [31 HI [5j [6] is that they all involve Hamiltonian 
constraints that satisfy an Abelian algebra. In the models in [H E] one is furthermore able to 
rewrite the Hamiltonian constraints in deparametrised form , that is in the form P + h = where 
P is the momentum conjugate to the time reference field and h is a function of all variables but 
the reference fields of density weight one, also called Kuchaf density. As a consequence it follows 
immediately that also the Kuchaf densities h mutually commute, that is {h(x), h(y)} = for all 
points x, y in the spatial hypersurface. For the reason that in the model in [6] the spatially 
diffeomorphism constraint is not rewritten in a form linear in momenta of reference field, the 
conclusion of strongly commuting Kuchaf densities h does not immediately follow but the algebra 
could close weakly, that is up to certain combinations of the spatial diffeomorphism constraint. 
However, as proven in [6] the corresponding Kuchar densities h of this model also commute 
strongly. 

Those Kuchaf densities are particular functions built from Q = y det(g), D = q ab C a Cb, C where 
q a b denotes the spatial three metric and C a , C respectively the non dust contributions to the 
spatial diffeomorphism and Hamiltonian constraints respectively (provided that matter couples 
only to the metric but not to its derivatives so that the hypersurface deformation algebra holds 
[30]). In general, rewriting the Hamiltonian constraints in this deparametrised form involves 
solving algebraic equations which leads to branches of the phase space labelled by the choice of 
certain signs, in this sense the description is only local on phase space and restricted to a single 
branch. 

The existence of the Kuchaf densities begs for the question whether one cannot obtain all 
functions h without going through an action principle. This problem was solved in [8j. Following 
[8] we introduce the density zero scalars d = D/Q 2 , c = C/Q and display a density w Kuchaf 
density in the form h = Q w K(c, d). Then the infinitely many equations {h(x), h(y)} = V x, y 
reduces, using the hypersurface deformation algebra, to the single first order PDE [8] 

A general integral can be obtained in terms of S = ln(K) by writing (|2.7p in the form 

wdS , /dS\ 2 1 fdS\ 2 2 . , 

which allows to separate the variables S(c,d) = Si(c) + 82(d) and to reduce (|2.8p to two simple 
quadratures which can be written in closed form and involve two arbitrary parameters. The 
complete integral can then be obtained by the envelope construction. 

All solutions that were obtained by considering models of type I or II are of course described 
by (|2.7p . One can ask the converse question how to build a Lagrangian of type I or II respectively 
which reproduces a solution of (]2.7|) . This was analysed in [9] and [10] respectively. In [9] it 
is shown that for all solutions of (|2.7j) a type I action involving a single scalar field T with 
non derivative coupling of the metric and a single Lagrange multiplier field p can be found 
reproducing it. Remarkably, all solutions of (|2.7p satisfying certain reality conditions (such as 
the positivity of the kinetic term) can be obtained from an action of the form 

£nrd = -\V\^i9)\ iP g» v U»U u + A(p)], U„ = V^T (2.9) 

That is, it is of the non rotating dust type just that the cosmological constant term is allowed 
to be a general function of the Lagrange multiplier. It has to satisfy an ODE that is matched 
to the two ODE's in (|2.8p . By contrast, the class of solutions of fl2.T|> that one obtains from the 
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type II models is more restricted. One finds from the Legendre transform (P = dL/dT) the 
following expression for p 2 := (P/Q) 2 as shown in [9] 

P 2 ■= (^) 2 = (L'(I)) 2 (I + q ab T a T b ) = (L'(I)) 2 (2-10) 

In the second step we used the spatial diffeomorphism constraint and applied the Brown-Kuchaf 
mechanism, that is using the relation q ab T^ a T^ = D/P 2 in order to replace the second term in 
the brackets above. The Hamiltonian constraint is given by 

ctot = Q { c+ W) ~ L{I) ) (2 - n) 

where / solves ()2.10p . To check whether a given solution p = —K(c,d) of (|2,7p corresponds 
to a Lagrangian L(I) we insert p = —K(c,d) into (|2.10p and solve the resulting equation for 
c = c(I, d). Then (|2.1ip must become the identity 

K(c(I,d),d) 2 

= c(I,d) + L(J) (2.12) 

which is an ODE for L provided the explicit dependence of (|2.12[) on d drops out. 



2.4 Global Parametrisation of Theory Space 

We now combine these known results and write down an action including up to 8 scalar fields 
and several parameters that describes all the models so far considered. Since the models of type 
II are contained in the set of models of type I in the fashion described we focus on the former 
set. Consider 

£d = --Vldet^W \P (V„T) (V„T) + A(p)V„V v + 2B(p)(V li T)V v ] + A(p)) , = W 3 V ^ 

(2-13) 

where A, B, A are arbitrary functions of the field p (we could also have considered an arbitrary 
function F(p) rather than p as the coefficient of the (VT) 2 term but this can be absorbed by 
a field redefinition F(p) = p' . Below we list the choices of these functions corresponding to the 
models reviewed in section (12. ip 

Ltd A = B = A = p 

C ND A = l, p = B = A = . . 

C NRD A = B = 0, A = p { ■ > 

C GD A = 0,B = l,A = p 

Among these models we find many new ones such as^4 = l, B = 0, A = p and of course all the 
models that one obtains by choosing A, B, A of a more general than linear form (notice however 
that unless A = we may absorb A into the Wj and have only B as a free function at our 
disposal. We keep both A, B so that we can treat the general case in a unified form). 



2.5 General Hamiltonian Analysis of Theory Space 

The general Hamiltonian analysis of the action corresponding to the sum of (|2,13p and the 
standard matter and geometry contributions proceeds roughly as follows (we consider the case 
that p and that not both of A, B vanish, more singular cases can be treated similarly, see 
[TT] for all the details in the case A, B ^ 0): 

Computing the momenta P, Pj conjugate to T, S 3 , the momenta it, ir 3 to p, Wj and the momenta 
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II, II a to lapse and shift functions N, N a in the 3+1 split of the action, we discover the primary 
constraints z = n = 0, z 3 = tt 3 = 0, Z = II = 0, Z a = as well as another set of 3 (2) linearly 
independent constraints Q 3 = for A ^ (A = 0) which involves only P,Pj,Wj (Pj,Wj), in 
both cases demand that PjWk — WkPj = and which come from the fact that one cannot 
solve for 3 (2) of the velocities S 3 . Altogether the canonical Hamiltonian depends on 11 (10) 
undetermined velocities. The stability analysis of the primary constraints with respect to the 
canonical Hamiltonian H C3uQ yields that one can solve 

{TT 3 ,H can }=0, j = l,..,3(2) {C J ,#can} = 0, j = 1, .., 3 (2) (2.15) 

for 6 (4) of the velocities and that there are the 5 (6) secondary constraints 

s = ^ = 0, C tot = ^ = 0, CT = = 0, (K = W", dH can } = 0) (2.16) 

The secondary constraints can be stabilised by solving for 1 (2) of the remaining unfixed velocities 
so that altogether 7 of 11 (6 of 10) have been fixed, leaving in both cased N, N a undetermined. 
Altogether 5 (6) secondary constraints were found leading to altogether 16 constraints. 

One finds that the 4 pairs (z, s), (z 3 , ( 3 ); j = 1,..,3 ((z, s), (z 3 , ( 3 ), j = l,2;(z 3 ,K)) form 
a second class system while linear combinations of Z, Z a , C tot , C* ot with the second class con- 
straints forms a first class set. The Dirac bracket between functions independent of p, Wj,N, N a 
and their conjugate momenta reduces to the original Poisson bracket. This follows from the fact 
that such functions have vanishing Poisson brackets with the constraints z, z 3 which mutually 
commute among themselves. The inverse of the Dirac matrix between the second class con- 
straints is therefore such that the difference between the difference between the Poisson bracket 
and the Dirac bracket involves at least one Poisson bracket with one of the z, z 3 . See section |4] or 
|17j for more details. One then solves the second class constraints strongly thereby eliminating 
p, Wj from the canonical Hamiltonian altogether. The result is the same as if one had eliminated 
right from the beginning p, Wj by using dH can /dp = dH C£LQ /dWj = and dropping the primary 
constraints from the Hamiltonian. This follows because these four equations are just four of the 
second class constraints. 

The canonical Hamiltonian is then a linear combination of the first class constraints C tot , C* ot 
which besides geometry and standard matter degrees of freedom, which we denote collectively 
by (q,p), now only depend on (T, P), (Pj,S 3 ) if not both A, B vanish and otherwise only on 
(T,P). Crucially, since in the first case the Hamiltonian depended only on the combination 
WjS J a and since Wj oc R on the constraint surface of the second class constraints, the canoni- 
cal Hamiltonian eventually depends only on the combination PjS 3 a as far as the S 3 dependent 
terms are concerned. This combination, however, is weakly equal to — (C a + PT a ) (see below) 
where C a is the contribution of geometry and standard matter to the spatial diffeomorphism 
constraint. It is for that reason (extended Brown - Kuchaf mechanism) that this particular 
combination was chosen in the original Lagrangian. In both cases, the explicit dependence on T 
in the Hamiltonian constraint is through T a only and the mechanism just displayed makes sure 
that this remains true also when substituting PjS 3 a - 

The remaining analysis of the system then proceeds as follows. It follows from the above 
that in the first case the constraints can be written in the equivalent form 

C tot = P + h(T,q,p), Cl ot = PT a + P,S 3 a + C a {q,p) (2.17) 

while in the second one has 

C tot =P + h(T,q,p), Cl ot =PT a + C a (q,p) (2.18) 
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In the timelike dust model the dependence of h on T even drops out completely while in (|2,18[) 
the fact that h depends only on q ab T a T b » q ab C a C b /P 2 can be used to further massage the 
constraints into the form 

C tot = P + %,p), C* ot = PT a + C a (q,p) (2.19) 

In all other models, h keeps an explicit dependence on T a . Nevertheless, in all models the 
Hamiltonian constraints in (|2.17p . (|2.18|) and (|2.19p strongly Poisson commute. This is because 
they are first class by construction but since the Poisson bracket eliminates the dependence on 
P, the bracket must actually vanish identically in the four reference field case (p ^ and not 
both A, B vanish). Notice, however, that this only implies that also the h strongly Poisson 
commute if they are independent of T. In the single reference field it is more complicated to 
show that the h strongly commute as has been discussed in |10] and [6]. 

In the case of model (|2.17p one can now perform a symplectic reduction of the spatial diffeomor- 
phism constraint by pulling back all tensors and spinors by the diffeomorphism x i— > a 3 = S J (x) 
which is a canonical transformation [lj [IT] to the effect that all degrees of freedom except Pj , S 3 
remain canonical pairs when expressed in the new frame while in the new frame the momen- 
tum conjugate to S 3 becomes ([dS ^/9x] _1 )"C* ot so that this canonical pair drops out from our 
attention. The remaining degrees of freedom are spatially diffeomorphism invariant when ex- 
pressed in this frame. We will denote them by (q' ,p' ,T' , P') in order to distinguish them from 
(q,p,T, P). In case of (j2. 19j) . spatial diffeomorphism invariant quantities must be constructed 
by other means. We will also pull back C tot (x) which then becomes the same function C tot '(a) 
of q',p'T',P' at a as C tot (x) was of q,p,T,P at x. 

The final reduction of the system now employs the projector formula discovered in I12| . 
See for a review in [T7] in the notation used here. Let F' , F be spatially diffeomorphism invariant 
functions on the phase space respectively which depend only on (q',p') in case of 112. 17|) and only 
on (q,p) in case of (|2.19p respectively. Let 

0' f ,{t) := (exp({C tot V], •}) • F%'=T'- T ; F (r) := (exp({C tot [g], .}) • F) 9=T . T ; (2.20) 

The notation means that we compute the Hamiltonian flow of the constraint indicated when 
smeared with a numerical test function g' or g and then set g' = T' — r or g = T — r. We cannot 
directly insert T — r into the exponent due to the phase space T dependence. Formula (I2.20|) 
is rather remarkable in several aspects. First of all (|2,20p strongly Poisson commutes with all 
constraints. Secondly, we have the equal time r Poisson brackets 

{0> fi (t),0' gi (t)} = 0> {f , g , y {t) (2.21) 

where {.,.}* is the Dirac bracket with respect to the second class pair T, C tot ' (similar statements 
hold for the unprimed quantities). For quantities F',G' independent of P' it reduces to the 
normal Poisson bracket. It follows that if F',G' are conjugate so are 0' f ,(t), 0' g ,(t). Secondly 
we have due to the mutual strong commutativity of the C tot/ (er) 



d_ 

7h 



°'f'^ = E(^T)j/ d '^-j d 3 a n [T-T\a 1 )Ur-T'(a n ^)]{C tot (a 1 ),{..{C tot (a n ),F'} 

= 0{c7t°t'[l],F'}( T ) = 0/i'[l],F'}( T ) = ^{h'[l],F'}*( T ) 

= {0' mi] (T),0' F/ (r)} (2 

where in the second step we realised the most inner Poisson bracket as the one between C tot/ [l] = 
/ d 3 aC tot '(a) and F' , in the third we used independence of F' on T' and that C tot/ = P' + ti, 
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in the fourth we used independence of both h', F' of P' and in the last we used (|2.2ip . The next 
property we need is that F' i— > 0' f ,(t) preserves the multiplicative and additive structure on the 
Abelian algebra of phase space functions, that is, 

0' f ,{t) = F\0' q ,(T),0' p ,(T),0' T ,(T)) (2.23) 

if F' = F'(q' ,p' ,T') depends on q',p',T' only. Now, the observables associated to the time 
reference field is given by 0' t ,(t) = r where r is a spatial constant. Likewise we obtain for the 
spatial derivative of T with respect to the (dust) spatial coordinates T a j the following observable 
0' T i CT (r) = d aJ T = 0. Accordingly, since h! only depends on the derivatives of T' because we 
assumed that we do not consider potential terms of the reference matter fields in all models and 
not on T' itself we can define the intermediate Hamiltonian by 

H> := J d 3 ah'(T',q',p') T , =0 (2.24) 

Then (pT22D becomes 

^0' f ,{t) = {0' h ,(t),0' f ,(t)} = 0' {h , jF , y {t) = 0> {h , }F , } (t) (2.25) 

where independence of H' of P' was used again. The properties of F' that were used in this 
derivation is that F' is independent of both T' , P' . Formula (I2.25P therefore in particular applies 
to F' = H' whence 

±O' h ,(t) = (2.26) 
is in fact r independent. It follows that (|2.25p can be rewritten as 

±O f ,(t) = {H,O f ,(t)} (2.27) 

where we have identified the physical, fully gauge invariant Hamiltonian 

H:=O' H ,(0)=O' m (0) (2.28) 
and in the second step we have made again use of O^,(0) = 0. 

Alternatively, one can derive the fact that the physical Hamiltonian H is not depending on 
time in the models considered as follows: 

In case the function F' depends explicitly on r the formula in (2.22) needs to be modified and 
we obtain 

T 0' f ,(t) = {0' hl[1] (r),0' F ,(r)} + (2.29) 
Now, let us consider the case that h'[l] depends on q',p',T' but not on P', then we have 

h'[l] = J d 3 ah'(T',q',p') (2.30) 

Using that the application of the projector F' \— )■ O f ,(t) preserves the multiplicative and additive 
structure on the Abelian algebra of phase space functions, we obtain 



H = 0' h[1] (r) = J d 3 ah'(0' T ,,O q ,,O p ,) (2.31) 



Assuming that we do not consider potential terms of the reference fields, as it is done in all 
models considered in this paper, they occur only with spatial derivatives. Now the observable 
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associated to the time reference field is just given by 0' t ,(t) = r. Since r does not depend on 
the spatial (dust) coordinates we further get 0' T , s (r) = r a j = 0. Consequently, in this case 
the physical Hamiltonian is of the form 

0' m {r) = J d 3 ah'(O q/ (r),O p/ (r)) (2.32) 

and has no explicit r-dependence. Knowing this, we can immediately see that it is constant 
with respect to the r-evolution using ()2.29|) 

^0' h/[1] (r) = {0' m {T),0' m {T)} + = (2.33) 

Therefore we can use any value of r for the fully gauge invariant Hamiltonian H and define 

H := O' h/[1] (0) (2.34) 



Using (|2.23p these are precisely the Hamiltonian equations of motion with respect to the Hamil- 
tonian (I2.34p . Notice that unless h/< T ,_ is a Kuchaf density, H' itself is not invariant under the 
Hamiltonian flow of generated by C tot/ , whence in this case the projector on H' in (12.24p is not 
trivial and H' must be replaced by (|2.28]h Also, unless h', T , =0 is a Kuchaf density, it will not be 

preserved by H whence H in general has only the three sets of local Noether charges Cj which 
arises as the pull back of C a by the diffeomorphism S J (notice that h', C totl are a scalar densities 
of weight one in a space and that r — T is a scalar so that (|2.20|) maps the diffeomorphism 
invariant quantity h'[l] to a diffeomorphism invariant quantity). 

We now reformulate in manifestly gauge invariant form. Consider the Dirac observables 
q = O' q ,(0), p = Op,(0), H = 0^(0) = h'[l](q,p) = H'(q,p), where we assume that h[l] 
includes only derivatives of the time reference field. Then (|2.2ip assures that Poisson brackets 
between functions of the hatted quantities can be computed by computing the brackets between 
the same functions of the primed quantities and then evaluated at the hatted quantities. 



3 Non Rotational Dust 

We carry out some of the details sketched in section (j2.5j) for the choice A = B = 0,A = p from 
a manifestly gauge invariant viewpoint and compare the results with [23]. Of course, almost 
everything we say here is already contained in [UET]- We will show that this model is contained 
as a subsector in the model |21| . both classically and in the quantum theory and also explain 
why proposal for the quantum theory in [24J cannot describe the physical sector of the model. 

In the case of non rotational dust the primary constraints are ir = U = U a = and the 
canonical Hamiltonian is given by 

#can = N a C™ + N \C + \ ( + p^/d^)[q ab T a T b + 1] ) ) +V7T + U U + U a U a (3.1) 

V 2 \PV det (i) J J 

with undetermined velocities v,U,U a . Stability of the primary constrains requires 
cT = PT a + C a = 

C tot = C + ~ ( + p^(q)[ q ab T a T b + 1] ] = 

P '~ + ^det(q)[q ab T a T b + 1] = (3.2) 



p 2 v / det(g) 
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Since {tt, s} ^ 0, the constraints S, tt form a second class pair and stability of s can be ensured 
by choosing v. We add to C* ot the term 

^ 'P,a proportional to the second class constraint tt. Then 
this enlarged constraint generates spatial diffeomorphisms on all variables contained in C^ * , C tot 
and thus preserves these secondary constraints. The smeared C tot Poisson commute to a smeared 
spatial diffeomorphism constraint according to the hypersurface deformation algebra because, 
up to the factors of p, the dust contribution C D to the Hamiltonian constraints coincides with 
that of a massless Klein-Gordon field plus a cosmological constant term, {C,C D } does not 
involve Poisson brackets between derivatives of fields and C D , C D lets its derivatives act on 
the smearing fields only so that the factors of p,l/p in C D cancel for the same reason as 
y / det(g), l/-\/det(g). It follows that there is only one second class pair. The corresponding 
Dirac bracket {/, g}* differs from the Poisson bracket by terms that involves terms of the form 
{tt, f}{n, 9}, {s, /}{vr, g}, {tt, f}{s, g} since {vr, vr} = / {s, s}, {s, tt} (the inverse of the Dirac 
matrix Ajj = {si,Sj}, s± = s, S2 = vr has the zeroes in complementary places as compared to 
the Dirac matrix itself). Hence, the Dirac bracket of phase space functions independent of tt,p 
coincides with the Poisson bracket. We solve the second class constraints strongly and obtain 

* = °' ^ = det( 9 )[l + g ^T a T b ] (3 - 3) 

Choosing a sign e for p we find that the canonical Hamiltonian is a linear combination of first 
class constraints 

Cl ot = PT a + C a = 

This can be massaged into 



C tot = C + e\P\Jl + q ab T a T b (3.4) 



C tot = C + e^P 2 + q ab C a C b (3.5) 



using the first equation in f|3.4[) . On the constraint surface, C has the sign — e (unless all 
quantities vanish). Inverting for \P\ we find 

|p|2 = C 2_ (3 6) 

which is constrained to be non negative. Thus, choosing a sign 5 for P we may replace f)3.5[) by 



c tot = p _ 6y j C 2 _ q ab CaCb = P -Sh (3.7) 

In the limit q ab C a C b -)■ (JES) and (03} become C + e|P| = P - 5\C\ = which is consistent for 
either choice of sign e, 5. It is customary to choose e = 1 = —5 (this is one of the four possible 
subsectors of phase space that one must pick in these local considerations). As H no longer 
depends on T it must be a Kuchaf density which maybe checked explicitly. We stress here that 
that in this model, as one can see explicitly from (|3.4p or (|3.5p . the sign of C is constrained, 
specifically — e is the sign of C. In particular, for e = 1 (positive kinetic energy term), C is 
constrained to be negative. 

The physical Hamiltonian becomes by the general considerations of the previous section the 
manifestly positive expression (whatever choice of e, 5 would have been made, (|3.8p has positive 
sign) 

H= f d 3 xh(x) (3.8) 



generating time evolution of functions of the form 

O f (t) = [exp({C tot [ 5 ], .}) • F] g=T . T (3.9) 
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where F is spatially diffeomorphism invariant and independent of T, P. The functions (|3.8p and 
(|3.9p are full Dirac observables of the theory and their algebra may be computed on the constraint 
surface C tot = C tot = 0. Since they are gauge invariant, their algebra may be evaluated on the 
gauge cut T = tq for some fixed tq = const., say tq = within that constraint surface. On that 
gauge cut we have C* ot = C a = and h = \C\ as well as Of{t) = exp([r — ro]{/t[l], .}) • F, in 
particular Of{tq) = F. The generator of time evolution is evidently still given by H = h[l] in 
(|3.8p but on the gauge cut its Hamiltonian vector field reduced to that of |C|[1] as follows from 

{H,f} = J d 3 x^(2C{CJ}-2q ab C b {C a ,f}-C a C b {q ab ,f}) Ca =o 

= j d 3 xsgn(C){C,f} = {\C\[l],f} (3.10) 

A different way to see this is to start from the formulation (|3.4p and to write the Hamiltonian 
constraint in the equivalent form 

r tot = p+ l C l - =-p + h (3.11) 

VI + q ab T a T b 

where the sign choice e = 1 = —5 has been adopted. The constraints ([3. lip are still Abelian and 
the machinery sketched in the previous section applies but h is no Kuchaf density any longer 
and h[l] is no longer a Dirac observable. However, the theory sketched there shows that the 
physical Hamiltonian is given by (notice that 0.(t) is now defined in terms of h while O.(r) is 
defined in terms of h) 

H = O~ h[1] (0) = O mi] (0) (3.12) 

since Or a (0) = 0. As the algebra of the F and \C\ is isomorphic to the algebra of the O_f(0) 
and 0\c\ (0) we arrive at the same picture at this fully gauge invariant level as in its gauge fixed 
version. Yet a third way to see this is to start from (|3,8p and to make use of the fact that h is a 
Kuchaf density. Then trivially H = O#(0) = H. On the other hand in O#(0) we may replace C a 
appearing in h = ^C 2 - q ab C a C b by C* ot as O Ta (0) = 0. Writing C = O c (0), C* ot = O C tot(0) 
etc. we have 

H = j d 3 x ^J[C] 2 -g^qtotqtot (313) 

and now a calculation similar to (|3,10p reveals that on the constraint surface we may replace H 
by Oi£7|h](0). We conclude that on the constraint surface all three approaches are equivalent. 

Notice that without the additional fields S"- 7 spatially diffeomorphism functions F are difficult 
to construct: Typically one will form them from scalar densities of weight one, built from the 
the fields q,p and their spatial covariant derivatives. These are no longer elementary functions 
such as q,p in terms of which a quantisation is typically easy. The equal time Poisson algebra 
of the observables Of(t) is isomorphic to that of the F since they do not depend on the P,T. 
However, since it is unknown how to construct a complete and independent set of such F and 
thus how to write their algebra in closed form, it is not possible to use the Of(0) or equiva- 
lently the F as a platform for quantisation and thus a full reduced phase space quantisation 
is impossible for practical reasons. One thus has to resort to a hybrid scheme such as in [22j 
where one starts from the phase space unreduced with respect to the spatial diffeomorphism 
constraint and quantises the O/(0) whose algebra is isomorphic!! to that of the /. The spatial 
diffeomorphism constraint is then solved in the quantum theory following |31j . However, while 

3 In [22] actually the / were quantised and the projector formula O. (r) was implemented in the quantum theory. 
The results are equivalent. 
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this equips the theory with a physical Hamiltonian H (for the matter model considered in [51 [6] 
rather than that model of non rotational dust considered as a special case of [1] used in this 
section), the task of implementing the Op(0) remains to be performed and the reservations of 
[23] are sustained. 

Coming back to the model of this section that was studied in |24j the spatial diffeomorphism 
constraint needs to be solved at the quantum level for the reason that likewise to |22] only one 
(time) reference field is involved in that model. Hence in order to describe the physical sector 
of the quantum theory, the same procedure as in [22] must be adopted as far as the quantum 
theory is concerned and the physical Hamiltonian needs to be quantised on the diffeomorphism 
invariant Hilbert space. However, the authors of [23] suggest a different procedure. Rather than 
quantising the theory and the physical Hamiltonian |C| on the diffeomorphism invariant Hilbert 
%diff space constructed in [31] and employed in [22] , they suggest to use the algebraic quantum 
gravity framework of [21] and to quantise q,p, \ C\ directly on the corresponding Hilbert space 
Haqg- Since that Hilbert space is based on abstract rather than embedded graphs, their view- 
point seems to be that therefore the spatial diffeomorphism constraint (and by the above gauge 
fixing viewpoint also the Hamiltonian constraint) can be considered as solved so that Haqg is 
actually the physical Hilbert space with the simplest possible physical Hamiltonian C on it. We 
make five comments about this: 

1. First of all, since really the physical Hamiltonian is \C\ rather than C there is secretly still 
a square root \C\ := VC 2 involved which makes \C\ no better than h = y 'C 2 — q ab C a Cb 
when quantising the theory. However, note that there is no principal problem with the 
square root because one can use the methods of [26] to handle it. The authors of [23] 
could of course take the point of view that they fix the sign of C classically and then 
quantise the corresponding part of the phase space in order to avoid the sign function 
at the quantum level. However, this is almost impossible to control in practical terms, 
since it requires detailed knowledge about the spectrum of C. The only practical way of 
ensuring the sign is to introduce the absolute value (that is a square root). Alternatively, 
the authors of [23] could argue to relax the typical energy conditions of classical General 
Relativity which however is not unproblematic. In any case, this discussion shows, that 
the claim that the physical Hamiltonian simplifies in the maximal sense and just reduces 
to C deserves further discussion By contrast, as we will discuss in the next section, there 
exists already a model [3], the Gaussian dust model, in which one can obtain physical 
Hamiltonian densities that simplify in the maximal sense and are just given by C without 
the need to fix signs at the classical level. 

2. As noticed previously in the literature (see e.g. |27] and references therein) the Hilbert 
space %aqg cannot be identified with the physical Hilbert space for models where the 
spatial diffeomorphism constraint is not solved already classically. We repeat the argument 
here for completeness: 

The physical interpretation of spin network functions in %aqg whose abstract graphs 
are finite subgraphs of an infinite abstract graph is obtained by embedding them into 
the given spatial manifold a under consideration. Any such embedding Y establishes an 
isomorphism, at least when a is compact, with a subspace of the kinematical Hilbert space 
TihQG °f LQG. In that Hilbert space, spin network functions based on diffeomorphic graphs 
are gauge equivalent. Hence must be their preimages in Haqg under Y. For instance, if 
the infinite algebraic graph is cubic then all Wilson loop functions around plaquette loops 
should be identified. 

The physical sector of the quantum theory in AQG is determined by the kernel of the 
extended master constraint, which includes, for the reasons mentioned above, also a master 
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constraint version of the spatial diffeomorphism constraint at the algebraic level. The 
quantisation of the (not necessarily algebraic graph-preserving) extended master constraint 
is presented in [26] using techniques from [25]. (Note that the continuum version of the 
spatial diffeomorphisms constraint used at the embedded level has no counterpart at the 
algebraic level). Therefore, when solving the Hamiltonian constraint classically, as it is 
done in [24J, one still needs to find the solutions of the operator corresponding to the 
classical expression 

M := [ ^n qa bCaCb (3.14) 

J VMq) 

in order to solve the spatial diffeomorphism constraint at the quantum level and thus 
describe the physical sector of the theory. The quantisation of such an operator was 
presented in [25, 26J, however in the case of the model in [24] this quantisation cannot be 
copied. The reason is the following: Assume that one is able to find the solution space 
of the (corresponding) operator in A3. 14j) . Then this solutions space is not left invariant 
under the action of the physical Hamiltonian because q ab C a Cb will not commute with the 
physical Hamiltonian densities given by C (or |C| respectively when the sign is not fixed 
at the classical level). Therefore one needs to find a quantisation of M such that the 
corresponding operator commutes with C modulo M. Such a quantisation^ so far does 
not exist in the literature and also has not been worked out in [2T 



3. The authors of |24J sketch the prospects that their model may have. The list of correspond- 
ing items is almost identical to that of [21] which predates [23]. Hence let us compare the 
models of [24] and [21] . The model of [21] is based on the general model in pQ rather than 
the special case of non rotational dust also discussed in [lj. It thus performs, in contrast to 
|24] , a complete reduction of the physical phase space since it contains also the scalar fields 
S 3 . The authors of [21] outline a quantisation of the algebra of full quantum observables 
and the physical Hamiltonian H, which is isomorphic to the simple algebra of the q,p, on 
both Hilbert spaces Hlqg and Haqg which rightfully can be called the physical Hilbert 
space in contrast to the model [23]. The price to pay is that the physical Hamiltonian is 
now given by the quantisation of the term 



H = J d 3 x ^\C 2 - q ab C a C b \ (3.15) 

The additional absolute value respects the fact that the quantity under the square root is 
constrained to be non negative in the classical theory and in this form (|3.15p allows for a 
well defined quantisation as discussed in |21j . Since the spatial diffeomorphism constraint 
reads 

= PT a + P 3 S 3 a + C a (3.16) 

even in the gauge T = r and S 3 = a 3 (|3.15p does not reduce to |C|[1]. This is because in 
this gauge (]3. 16|) becomes 

CT = m + C a (3.17) 

However, this is neither a problem, since semiclassical tools are available in order to deal 
with square root Hamiltonians [26], nor worse than working with |C|. 

4. It is easy to see that the model [24] corresponds to a subsector of the model [21]. This can 
be demonstrated both at the classical and the quantum level. 



4 Note that it is not enough to commute modulo C a which classically is trivially the case because C a in contrast 
to M is not a well defined operator. 
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At the classical level this corresponds to the observation that the portion of phase space 
where C a = is left invariant by f)3. 15|) . Let us discuss the degrees of freedoms in both 
models. In [24] one has gravity plus one scalar field and in [21] gravity plus four scalar 
fields. In both models the total spatial diffeomorphism constraint is required to vanish 
yielding at the level of the gauge T = r for the model [24] the constraint C a = 0. When 
we consider for the model in [21] the gauge T = t and = a 3 we obtain C a + PjS J a = 
where Pj are the momenta conjugate to S 1 . Let us define P a := PjSi then in the model 
|21] we have C a = —P a and it can be shown that P a is a constant of motion and thus so 
is C a . In particular, requiring C a to vanish in the model of [21] means P a = yielding 3 
additional constraints reducing the number of degrees of freedom to that of the model [24J. 
Further evidence is obtained by noting that P a = means Pj = and from this follows, 
using Wj ~ Pj, that the Lagrange multiplier fields Wj are constrained to vanish. This 
means that the reference fields S 3 are no longer present in the action. Hence we arrive 
at the special case of non-rotating dust, demonstrating the embedding at the classical level. 

At the quantum level, the embedding is less intuitive but nevertheless can be shown. Let 
us first discuss the case where we choose the version of [21] where one quantises on %lqg- 
One observes that the physical Hamiltonian operator H is invariant under diffeomorphism 
on the dust manifold S(a) and thus commutes with it generator C a where the index a here 
needs to be understood as labelling coordinates on the dust manifold S(a). One would 
rather call C a a generator of symmetries in dust space than a generator of gauge trans- 
formations. This symmetry generator is implemented on the physical Hilbert space which 
is the standard kinematical Hilbert space %lqg of LQG. Now given the physical Hilbert 
space, one can look for functions that are invariant under the symmetry transformations 
generated by C a and those (dust) diffeomorphism invariant 'functions', that are rather 
distributions, will be the 'subset' of functions that are additionally annihilated by C a and 
hence also by the operator q ab C a Cb when q ab C a Cb is properly quantised as to annihilate 
diffeomorphism invariant distributions. On those diffeomorphism invariant distributions 
H reduces to |C|. 

Now, let us consider the second possibility where one chooses the version of |21] that is 
quantised on Haqg and see how the quantum theory of [2J] is embedded in that of [2T] . 
Suppose we have fully reduced the Hamiltonian and spatial diffeomorphism constraint 
at the classical level, then the operator associated with the classical expression shown 
in (13.14p is a generator of symmetries rather than gauge transformations at the physical 
Hilbert space %aqg- Likewise to the standard LQG quantisation, one can now look for 
'functions' (distributions) that are annihilated by that operator and consider the solution 
space as a 'subspace' of the physical sector of the theory. However, here in general this 
'subspace' will not be left invariant under the dynamics unless, as discussed before the 
operators M, C (|3. 14[) are quantised in such a way that they mutually commute. This 
would in any case be a desirable feature of the theory but this is not granted by the naive 
quantisations presented in |26] and employed in |24j . In other words, whenever M is de- 
fined in a satisfactory way in the model [24J at the algebraic level then it corresponds to 
a subsector to the algebraic version of [2T] . 

Finally, we observe the following curiosity when working in the standard LQG framework. 
Let us consider the 'subspace' of functions that are invariant under the symmetry transfor- 
mation generated by the C a on Hlqg- Now this 'subspace' has mathematically the same 
structure as the standard diffeomorphism invariant Hilbert space of LQG. Therefore, on 
this 'subspace' the physical Hamiltonian reduces to an expression where all terms involv- 
ing q ab C a Cb can be neglected and furthermore a graph-modifying quantisation could be 
applied. Hence, in this sense, in this subsector knot class modifying models can be redis- 
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covered although it would be a bit artificial to force a knot class modifying quantisation 
in that particular 'subspace' whereas on the 'rest' of the Hilbert space a graph preserving 
quantisation is adopted. 

5. A final comment concerns the claim made in |24] that their physical Hamiltonian is quan- 
tised free of "anomalies" as compared to [26]. Just to avoid confusion, let us try to interpret 
this statement. If what is meant that one single Hamiltonian operator (as compared to 
an infinite number of Hamiltonian constraint operators) commutes with itself, the state- 
ment is empty. If what is meant is that the Hamiltonian densities (as compared to the 
Hamiltonian which is the integral of the densities) mutually commute which is equivalent 
to the mutual commuting of the Hamiltonian constraint operators, the statement is wrong 
at the algebraic level in which the authors of [24] are working. The commutator of two 
Hamiltonian densities is not an algebraic version of the diffeomorphism constraint and in 
that sense there is an anomaly. This is again the reason why in [26J a master constraint 
approach towards all constraints has been adopted. If what is meant is that within the 
embedded LQG framework the Hamiltonian constraints commute on diffeomorphism in- 
variant contributions, this is correct when using the quantisation of [32] but not when using 
the algebraic quantisation of [26] that the authors copy. In any case the quantisation of 
|32] is inappropriate for a physical Hamiltonian rather than an infinite number of Hamil- 
tonian constraints because the resulting Hamiltonian would not even be symmetric in that 
quantisation scheme (for constraints, symmetry is not necessary as one is only interested 
in the joint kernel). Rather, as pointed out in [21J, one must use the graph preserving 
quantisation scheme (and a symmetric ordering) in order to obtain a symmetric operator 
which is spatially diffeomorphism invariant. Finally, as already pointed out in |21j . the 
issue of anomalies is much less critical in this reduced phase space quantisation approach 
as compared to the operator constraint quantisation approach because the number of de- 
grees of freedom has been correctly reduced already at the classical level. In that sense 
there can be no anomaly. However, one might still be interested in implementing the clas- 
sical hypersurface deformation algebra as a physical principle to reduce the quantisation 
ambiguities and this remains true for all models considered so far. 

4 Gaussian Dust 

It transpires that an optimal model would be such that 1. the physical Hilbert space is the usual 
%lqg (or Haqg) and 2. the physical Hamiltonian density is equivalent to just C and not |C| 
or \J\C 2 — q ab C a Cb\. None of the models [211 1221 121] has both features. Remarkably, we find 
exactly such a model in [3]. In what follows we sketch the classical treatment of this model in 
some detail since it differs slightly from that of [17} [T8] . In contrast to [3] we will not simply 
eliminate the Lagrange multiplier fields by their equations of motion but rather go through a 
careful Dirac treatment as outlined in section 12.51 The quantisation of this model will be copied 
from |2H [32] and we will therefore be brief on that point. For the physical interpretation of 
this model as a Gaussian reference fluid we refer to [3]. We just mention here that variation of 
the action with respect to p, Wj respectively leads to the conditions U^U^ = — 1, U^Vj 1 = on 

the vectors = g^T tU , Vj := g^ h ' S] v and the equations of motion for T, S J yield the geodesic 
equation S/jjU = and the conservation equation V ^{WjU^) = 0. Hence the integral curves of 
U, labelled by 5 J (x) = aj = const, describe an observer in geodesic motion and the vectors Vj 
are orthogonal to the corresponding T = const, surfaces. Hence in the corresponding reference 
frame the metric assumes the Gaussian form with unit lapse and zero shift. The energy mo- 
mentum tensor turns out to be T^ v = pU^U u + ^ • WjVj^U u ^ which has no spatial trace with 
UfJJv + g^ u and thus no pressure whence the name "dust" is appropriate. 
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The Gaufiian dust action reads explicitly 

£ GD = -| det(g)\ l / 2 {p±[g^T,T u + 1] + g^TjW^]} (4.1) 

Performing the 3+1 split in the ADM frame we obtain 

gTT^T,,, = —[L n T] 2 + q ab T a T b , g^TjWjS^} = - [L n T] [WjL n S 3 ] + q^T^W^}; 

L n = n = ^[d t -N a d a ] (4.2) 

where N, N a are the usual lapse and shift functions and q ab is the intrinsic metric on the 3- 
manifold a (with inverse q ab ) which is mapped via a one parameter family of embeddings Yj into 
into a one parameter family of spacelike hypersurfaces = Yt{a) that foliate M. The timelike 
vector field n is unit normal to the foliation and its action on the scalars in (|4.2p coincides with 
the Lie derivative L n . 

Performing the Legendre transform we find together with -^/| det(<?)| = Ny/ det(g) 



C 



P := -^ = ^fettij{p[L n T] + W j [L n &}} 
oT 



C 



Pi ■= -^j = VMq)W j [L n T] 

^GD n 

op 



(4.3) 



The detailed constraint analysis can be found in the appendix. From it follows that we have 
8 first class constraints Z, Z a , C tot , C* ot and 8 second class constraints z,Zj,Ci,s,K. Fortu- 
nately it is not necessary to compute the corresponding Dirac bracket explicitly by the following 
argument: 

We arrange the second class constraints into 2 sets {K^\x)}j =l = {z, z-*}(x) and {K^\x)}j =l = 
{Ci, s, K}(x). Then the Dirac matrix 

A I l / u (x,y):={K^(x),Kj / J \y)} (4.4) 

and its inverse assumes a square block structure of the form 

A n = A 12 \ ! _ / ([A 12 ]" 1 ) T A 22 [A 12 ]- 1 -([A 12 ]- 1 ^ \ 

_[ A 12]T A 22 J ^ " ^ [A 12]-1 J 

As we will eventually solve the second class constraints for the Lagrange multiplier fields 
7T, tt 3 , p, Wj we are interested only in the restriction of the Dirac bracket to functions / of 
the variables q,p;T, P; , Pj. For such functions we have {K^x),/} = 0. Since the difference 
between the Dirac bracket and the Poisson bracket between /, f contains, due to (|4.5j) . only 
terms with at least one of {K^x), /} or {K^x), /'}, the Dirac bracket and the Poisson bracket 
actually coincide on the functions of interest. 

The solution of the second class constraints (j = K = s = is given by 

Wi = W 3 ^ 



w 3 



p 3 

p2 

2 ^3 



Q 2 (l + q ab T a T b ) 

(4.6) 




q^TaPS] 



1 + q ab T a T b 
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The solutions for Wj, p are not explicitly needed but we need to choose a sign e for W3/P3, insert 
the root of the second equation in f)4.6j) into the Hamiltonian constraint and see that the term 
proportional to p vanishes identically (since it enters linearly into the Hamiltonian constraint) 

<i'' lr i:,,i',s ! . 



C-tot = C + e[ PJl + q ab T a T b + -. ''"7^ ) (4.7) 
V V V 1 + Q ab T a T b ) 

while the spatial diffeomorphism constraint becomes 

C£ ot = PT a + PjS% + C a (4.8) 
Substituting for PjS 3 A in (|4.7p we find the equivalent and simple form 

C^ = C + e^= ~ qabT ' aCb (4.9) 
VI + q ab T a T b 

or equivalently 

ftot D U . D 1 /l 1 T' „l6 



c tot = p + ^ = . p + eC ^ 1 + ^oby^ _ ^"T a a (4.10) 

We see that h is of the type described in section (|2.5p so that the general theory outlined there 
applies. Moreover, we see a crucial difference with the model described in section [3) The sign e 
of W3/-P3 is unrelated to the sign of P, C. The choice of e in the Lagrangian also has no physical 
significance since the Wj dependent term is neither bounded from above nor from below. For 
definiteness we simply choose the phase space such that e = +1. Hence there is no absolute value 
of C involved in (|4.10p and the physical Hamiltonian becomes simply, following the exposition 
of section 12.51 

H = j d 3 a 6(a), C{a) = C(q,p) (4.11) 

where S = S(a) is the dust particle manifold and q = O' ql (0), q = 0^,(0), that is, it is the 
same function of q,p as is C of q,p. Here q',p' denotes the pull back of the fields q,p by the 
diffeomorphism a J = S J (x) and O'(0) is the projector map defined in section [231 The Poisson 
algebras of the q,p and of the q',p' are identical. Notice that (|4.1ip is invariant under Diff(5). 

The quantisation of the system can now be copied from |32[ 121] , If we work on the Hilbert 
space %lqg an d want the quantum operator to have the symmetries of its classical counterpart 
then H must preserve the subspaces %lqg,7 defined by the closed linear span of SU(2) invariant 
spin network functions over the graphs 7 embedded in S in order to be densely defined. In [21] 
it is described how this can be achieved using the notion of a minimal loop attachment and a 
corresponding projection operator onto "Hlqg,7- * ne framework of %aqg the operator H does 
not need to preserve any of the subgraphs of the algebraic graph and can be defined in terms of 
next neighbour loops |21l [26] just as in lattice gauge theory. In both approaches a symmetric 
ordering of the Hamiltonian densities must be chosen. Here the following subtlety arises: 
When considering \C\ rather than C the ordering of C is chosen in such a way that C acts only 
at the vertices of a graph. This ordering is not symmetric but this does not matter because one 
considers the positive operator valued density \C\ := \J C*C. However, when quantising C itself, 
such an ordering is not available. To make sure that C is symmetric and is densely defined (the 
danger being that the curvature term involved in C acts everywhere, not only at the vertices of 
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a spin network function) one uses the tools developed in [32] and writes 



(det(e(x))) 2 = f det({A(x),V t (x)}^ 
^det(q(x)f V Vdet(?(x) , 

fdet(Tr(A e (x)- 1 {A e (x),V e (x)})' 



lim 



^ J 



(4.12) 



— | lim 

2 / e->0 

where e(x) is the cotriad, A{x) is the connection, V e (x) is the volume function of a neighbour- 
hood of x with coordinate volume e 3 and A e (x) is a set of three SU(2) valued holonomy functions 
along three edges starting at x with linearly independent tangents at x and which span a coordi- 
nate volume e 3 as well. As shown in [32], each of the two factors in (|4.11j) admits a well defined 
quantisation in the limit of vanishing regulator and acts only at vertices. We can therefore 
freely position one factor each from (|4.11j) to the outmost left and right of the operator and then 
order symmetrically (including the projections on Hlqc-j m the case of a quantisation on Hlqg)- 



We notice that since C is not a Kuchaf density, the physical Hamiltonian has less symme- 
tries than for the model described in [17} l2l] , it does not Poisson commute with its density and 
thus it only has the Noether densities Cj(a) as conserved charges. However, C(a) becomes a 
Noether charge when the other three Noether charges Cj(a) vanish because the C(a), Cj(a) 
obey the hypersurface deformation algebra. Thus, on this sector of the classical theory, the 
results of [UJ El E3] continue to hold which in the case of [T7J [18] only assumed C(a) to be 
preserved and non vanishing while Cj(a) was kept arbitrary and in the case of [33] one was 
working with Cj{a) = anyway. As a result, the model described in this section is in agreement 
with the usual gravitational waves, cosmology, cosmological perturbation theory and black holes 
description as as described by geodesic test observers, following the analysis of [T7 ] [TB I [21] . 



5 Summary and Outlook 

In this paper we have accomplished three tasks: 

1. We have described all matter models considered so far for the purpose of deparametrisa- 
tion of General Relativity with an eye towards Quantum Gravity starting from a general 
Lagrangian. This serves to bring order into the space of models studied already and those 
that can be within the theory space described. We have described the Hamiltonian analysis 
of these models in a uniform fashion. 

2. Basically two types of models have been constructed: II. Those which deparametrise only 
time and I. those which deparametrise spacetime. 

We have shown that those of type II. necessarily must implement spatial diffeomorphism 
invar iance in the quantum theory. This is a non trivial subject for several reasons: 
First, even though the unitary implementation of spatial diffeomorphisms in LQG is no 
problem, there is no unique Hilbert space structure on the diffeomorphism invariant distri- 
butions [31] . To fix the ambiguities, the implementation of the ""relations among spatially 
diffeomorphism invariant observables may be consulted as a guiding principle, however, 
neither is a complete and algebraically independent set of such objects known even at the 
classical level nor is it clear that these can be quantised as to yield an honest representa- 
tion of their * algebra (see [34] for the tremendous difficulties encountered already in the 
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much simpler setting of the closed bosonic string). This is due to the fact that spatially 
diffeomorphism invariant functions on phase space are typically not simple polynomials of 
the elementary fields. The considerations made in |23j indicate that such anomalies are a 
conceivable possibility and that no choice of the ambiguity parameters may exist for which 
the representation problem can be solved. 

These difficulties are naturally avoided in the type I models where also the spatial diffeo- 
morphism symmetry is reduced already at the classical level. The non trivial observation 
is that this can be done while keeping a very simple Poisson * algebra of basic fields so that 
one can easily find Hilbert space representations thereof. This also extends to the reduc- 
tion of the symmetries generated by the Hamiltonian constraint, at least for the matter 
models considered here. This is quite remarkable when recalling the difficulties that one 
meets when trying to implement the Hamiltonian operator constraint in LQG, see e.g. the 
discussion in [2H [261 ESI [35] . 

3. The authors of [23] point out that their model, which belongs to type II, avoids a square 
root in the physical Hamiltonian density in contrast to the models considered in |21l [22] . 
We have shown that this statement has to be made more precise. If one insists on the 
usuall energy conditions of classical General Relativity then the analysis shows that rather 
\C\ = yCP 1 is the physical Hamiltonian density which again needs a square root. We 
remark that square roots are a nuisance rather than a caveat as semiclassical tools are 
available that can deal with this problem [26J. However, it is certainly true that the analysis 
simplifies significantly when there is no square root. Unfortunately, in type II models the 
quantum reduction of the spatial diffeomorphism constraint still poses a challenge. In view 
of what was just said, this leads to the natural question whether there are matter models 
within type I which lead to a physical Hamiltonian without square root. It is a remarkable 
achievement of Kuchaf and Torre [3] to have achieved just that in the form of the Gaussian 
dust. We have applied essential steps of the analysis performed in [T7] to the model [3] in 
order to show that all that was said in [171 118[ [2Tj 133] continues to be correct. We believe 
that in the class of models of type II considered, the degree of simplicity of the physical 
Hamiltonian cannot be improved over the model presented in section [H 

We believe that switching from the earlier operator constraint reduction (Dirac) approach that 
dominated the research in LQG over the past 20 years to the reduced phase space approach, 
concretely implemented within LQG for the first time in |17[ I2T] has many merits. Within 
the former approach, there are many steps of significant technical complexity to be overcome 
before one reaches the same level of technical and conceptual control as in the latter. After 
all, ultimately also in the former approach one aims at the physical inner product, physical 
vector states, physical observables and the physical dynamics which come for free in the latter 
approach. Moreover, on physical grounds the intrinsic description of the physics in terms of 
material reference systems makes a lot of sense, hence the approach advocated here and by 
many others before is well motivated. 

There are certainly some possible caveats: For instance, the classical reduction assumes that 
certain objects such as det(<9Sy<9x) is everywhere non vanishing. As one can show, this is a gauge 
invariant condition but it imposes a non holonomic constraint on the phase space. Hence, the 
description is not valid in the full unconstrained phase space. On the other hand, these issues 
are maybe less important than one thinks as the material reference system serves a role quite 
similar to three of the four real degrees of freedom of the Higgs field which give mass to and 
Z bosons: They are simply absorbed into those vector bosons thus making them local Isospin 
SU(2) singlets. Similarly here, the dust field is absorbed into the geometry and matter degrees 
of freedom, thus making them local Diff(M) singlets. Hence the material reference system is the 
gravitational counterpart of the Higgs effect. 
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The real physical question is whether scalar fields such as the Higgs field or the dust con- 
sidered here exist in nature. The answer might quite well be negative. Although the dust 
field considered here and elsewhere in some sense make it an ideal dark matter candidate since 
it is only interacting with gravity and with itself, it may nevertheless be phenomenologically 
excluded. What we need is maybe a realistic dark matter candidate, to be found in some ex- 
tension of the standard model, that can serve the purpose of deparametrisation. It would be 
perhaps most economic to isolate four degrees of freedom from the spacetime metric tensor for 
the purposes of deparametrisation but this is more difficult than with matter because one needs 
to construct scalars using covariant derivatives which makes such a construction non local and 
thus impractical. 

In any case, dust fields and their generalisations as described in the present work provide a 
proof of principle that a reduced phase space quantisation approach to LQG not only works but 
moreover has many advantages over the Dirac operator constraint approach because it lifts all 
the structures previously found in LQG directly to the physical versus kinematical level. In view 
of these advantages, the other caveats mentioned are in our mind problems of lower priority. 
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A Constraint Analysis for the Gaussian Dust Model 

In this section the constraint analysis of the Gaussian dust model is presented. From the 
Legendre transform we deduce the following primary constraints 



Z = 7T, Z J = 7T J 



z = n, z a = n , w -.= j kl p k Wi (a.i) 



with the momenta II = n o = for lapse and shift functions respectively. We notice that only 
two of the 3 constraints R? are linearly independent because obviously WjBP = 0. In what 
follows we will assume that ■ Wj > and £^ • Pj > 0- Then in a patch with W3, P3 / we 
can choose 

Cl := Pi ~ tt^-Pz = ~R 2 /W 3 , C2 := P2 ~ t^Pz = R l /Wz (A.2) 

as the independent ones which tell that all Pj are proportional to P3. The velocity combinations 
that can be solved for in (14.31) are 



[L n T} = P , Wj[L n S j ] = , P P P3 (A.3) 
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We compute the canonical Hamiltonian density using the abbreviation Q = \J det(g) 

-^GD.can = [PT + Pj S j +7Tf> + Wj + UN + U a N a - C GD ] ^ 

= {N a [PTa + PjSi] +uz + Uj z j + VZ + V a Z a + v 1 ^ + N(P[L n T] + ^[WjL f 

W 3 

+NQ{ l -p[-[LnTf + q^Ti + 1] + [-[L a T][Wj[L n S']] + q" t T a \W,S' J }} {m 



N-lPTa 4- P.Sj, ] + «, + + VZ + VZ a + „' & + N( P J^ + * [| - ||]) 

D D „ p 

A^PTa + P^] + + ^ + VZ + F a Z a + + - ^[-^f] 

+NQ{±p[q ab T a T b + 1] + q ab T a [W 3 S j b ]} 

N a [PTa + P^J +uz + Uj z j + VZ + V a Z a + [v 1 - N^-Qq^TaS^d 

Pz 

l W 3 Q 2Q [ W 3 l ] 



+NQ{-p[q ab T a T b + 1] + -^ q ab Ta[PjS\}} (A.4) 

where we used the abbreviations 

u = p, Uj = Wj, V = N, V a = N a , v 1 = L n S J ; I = 1, 2 (A.5) 

The contributions to the canonical Hamiltonian from geometry and standard matter are N a C a + 
NC which has to be added to (|A.4j) in order to obtain the total canonical Hamiltonian 

#can= / d 3 x(U GD ^ n + NC + N a C a ) (A.6) 

Stability of the primary constraints yields on the constraint surface of the primary constraints 

i pp n P 1 W 

@a := {Za> #can} = C a + PT a + PjS*^ 

:= { z, J F/ can } = y{-i[^] 2 + Q[^T a T fe + l]} 

{*,,fW = ~W ~ N^Qq^TaS 1 ,]^- 

{C/,#can} =U I ^-+M I (A.7) 

where we used that {Cf>C/} = an d Mj is some complicated and non vanishing expression 
whose explicit form will be of no further interest. The set of equations (|A.7p has to vanish which 
is accomplished by choosing 

v 1 = N^Qq ab T a S* b , uj = -§-Mj (A.8) 
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and imposing the secondary constraints C tot , C* ot , s, K. Notice that due to (|A.8j) K simplifies 
to 

pp n p2 1 

K = + + Q—q ah T a [PjS 3 h } (A.9) 

and the terms proportional to (j drop from the Hamiltonian. 

Turning to the stabilisation of the secondary constraints, we add to C* ot the linear combina- 
tion of the already stabilised primary constraints 7rp iQ + ^ Wj, a + IliV a + U b N b a + {U a N b ) :b . Then 
Qtot generates spatial diffeomorphisms on all variables involved and since T^can is a spatial scalar 
density of weight one, C* ot is stabilised. Next we have on the constraint surface determined so 
far 

= {C tot , J F/ can } = {C tot ,C tot [iV]}-n^-u 3 



dp dW 3 
|,7tot )C ,tot^| _ ^ S _ U 3 K 

{C tot ,C tot [N]} 



= {s,H can } = -2u^ + M 

= {K,H CSbIl } = -u^^+M' (A.10) 

where M is independent of u and M' depends linearly on u s . We can therefore solve the two last 
equations for u, u 3 respectively so that s, K are stabilised. As far as the first term is concerned 
we write for some smearing function / 

C GD [/] := / d 3 xf(T + U) (A.ll) 
:= / d 3 xf{±(^-^] 2 ) + Q(±pl q ab T a T b + l} + ^ 

and similar for C[f]. Then 

{C tot [f],C tot [f']} = {C[f],C[f}} + {C[f],C GD [f'}} - {C[f'},C GB [f]} + {C G v[f},C G v[f}} 

(A.12) 

The first term gives —C a [q ab (ff' b — /&/') as is well known from the hypersurface deformation 
algebra [30]. The second and first term cancel each other because the only piece from C that 
contributes is the gravitational piece which acts ultralocally only on the q ab dependence in Cqd 
which thus gives due to the non derivative coupling a term proportional to 8(x,y)[f(x)f'(y) — 
f(y)f'{x)\ = 0. The last term is of a new type. Again due to vanishing ultralocal terms we 
just need to focus on terms that lead to derivatives of the 5 distributions. For the same reason 
we only need to keep track of the x, y dependence of the smearing fields. Accordingly in the 
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following calculation we neglect ultralocal contributions 

{C GD [f],C GD [f']} = I f d 3 x fd 3 y[f(x)f'(y)-f(y)f / (x)}{C GD (x),C GB (y)} 



d 3 x / d 3 y[f(x)f'(y) - f{y)f{x)\ 



2 
1 
2 

x [{T(x),T(y)} + {T(x),V(y)} - {T(y),V(x)} + {V(x), V(y)}] 
{T(x),T(y)} = =0 



{T{x),V(y)} = ^{(PP 3 )( x ),(ip[^r a r 6 + l] + ^rt a [P i ^])(y)} 

-^{P 3 (x), (^fe a6 r a T 6 + 1] + ^q^TaiPjS^y)} 
Pi , W-i P u 



= ^g ab [P,^ a + PT a ] 
{V(x),V(y)} = QV^q^TdWlKPjS^/Ps^^P^/P^y)} 

= QV b T b q cd T d WiS^(x), S k d (y)}^ -x^y 

^3 P3 

= Q 2 q ab T b q cd T d WiS; b &x), S k d {y)}^ -x^y 

P3 P3 



Q^T^TaWiS^mixl S^y)}-^ - {^(z),^)}-^] - s y 

*3 ^3 

(A.13) 



It follows 



{C GD [f],C GD [f}} = I d 3 x J d 3 y[f{x)f{y)- f(y)f(x)]{T(x),V{y)} 

d 3 x [ d 3 y[f(x)f'(y) - f(y)f'(x)]5^q ab C GBa 



= -C G D[q ab (ff: b -f, b f')} (A.14) 

and therefore 

{C tot [f}, C tot [f'}} = -CT[q ah {ff[ b - /, 6 /')] (A.15) 

satisfies the hypersurface deformation algebra. All constraints are now stabilised. 

The constraints Z, Z a are trivially first class. The constraint C a ot is first class as we have 
already seen due to its geometrical action. The constraints I = 1,2 form second class 

partners as well as (z, K) and (213, s). Finally, the constraint C tot closes to itself with Z, Z a , C* ot 
and with z,z 3 respectively it closes to s,K respectively. To make it close with zi,Ci,s,K as 
well we add suitable linear combinations of all second class constraints to C tot so that it has 
vanishing Poisson brackets with all second class constraints. 
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